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We address the issue of unbiased model selection and propose a methodology based on Bayesian 
inference to extract physical information from kaon photoproduction p(-y, K + )A data. We use 
the single-channel Regge-plus-resonance (RPR) framework for p(-y, K + )A to illustrate the proposed 
strategy. The Bayesian evidence Z is a quantitative measure for the model's fitness given the world's 
data. We present a numerical method for performing the multidimensional integrals in the expression 
for the Bayesian evidence. We use the p(pf, K + )A data with an invariant energy W > 2.6 GeV in 
order to constrain the background contributions in the RPR framework with Bayesian inference. 
Next, the resonance information is extracted from the analysis of differential cross sections, single 
and double polarization observables. This background and resonance content constitutes the basis 
of a model which is coined RPR-2011. It is shown that RPR-2011 yields a comprehensive account 
of the kaon photoproduction data and provides reasonable predictions for p(e, e'K + )A observables. 

PACS numbers: 14.20.Gk, 14.40.Df, 11.55.Jy 



I. INTRODUCTION 

How to extract the nucleon resonance (N*) con- 
tent of the open strangeness photoproduction reactions 
p("f, K + )A is a long-standing question. Various analyses 
lead to disparate outcomes concerning the set of reso- 
nances that are likely to contribute [THE]. The recent 
availability of abundant high-statistics data has not pro- 
foundly changed the situation so far. This indeterminacy 
for the open strangeness channel is in stark contrast to 
the situation for pionic channels, where the contributing 
resonances can be successfully identified by means of a 
partial wave analysis for invariant energies W < 1.8 GeV. 
In open strangeness channels, this technique is less pow- 
erful as the nonresonant, or background, contributions 
are larger. The importance of background contributions 
calls for a framework which accounts for resonant and 
nonresonant processes and which provides a means to 
constrain both classes of reaction mechanisms indepen- 
dently. 

An efficient way of pinpointing the background am- 
plitude involves Regge phenomenology [HJ QH] ■ We will 
describe the p(^, K + )A reaction in the so-called Regge- 
plus-resonance (RPR) model, which combines ingredients 
of Regge phenomenology with elements of a typical iso- 
bar approach. The latter belongs to the class of tree-level 
effective Lagrangian models. In the RPR framework, the 
background amplitude is constrained by optimizing the 
adjustable parameters of a Reggeized background model 
to data obtained at sufficiently high energies so that the 
contribution of individual resonances is projected to be- 
come marginal j9[ [11] . 

Even with a properly constrained background contri- 
bution, the identification of the contributing resonances 
to p(j, K + )A remains a precarious task. Adding res- 
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onances increases the amount of adjustable parameters 
and improves the quality of the fit to the data. It stands 
to reason that one should not add more resonances than 
strictly necessary, in order to obtain a good model. One 
of the guiding principles for model selection is Occam's 
razor 12 . This principle dictates that if one has to 
choose between a simple and a more complex model, all 
else being equal, the simpler one should be preferred. In 
a realistic situation, however, all else is not equal, and 
this guiding principle should somehow be translated to 
a quantitative measure which balances between model 
complexity on the one hand and accuracy on the other 
hand. Such a measure can be derived from first princi- 
ples using Bayesian inference. This measure, called the 
Bayesian evidence Z, evaluates the overall performance 
of the model while penalizing for excessive complexity. 

In recent years, much effort has been directed to- 
wards a more comprehensive description of both elec- 
tromagnetic and hadronic meson production reactions 
from the nucleon within coupled-channels frameworks 
Uni [H] . Ideally, one would like to apply Bayesian in- 
ference to a state-of-the-art dynamical coupled-channels 
model. However, due to the multidimensional integrals 
involved in the computation of the Bayesian evidence, 
one is stricken by the curse of dimensionality: at worst, 
the computational cost increases exponentially with the 
number of adjustable parameters. Even if the number of 
adjustable parameters is kept in check, the sheer num- 
ber of model evaluations required for a Monte Carlo in- 
tegration calls for a realistic model of modest complex- 
ity. In this work, we will apply Bayesian inference to 
the single-channel RPR framework for kaon photopro- 
duction. We will consider several variants of the RPR 
model and use Bayesian inference to select the most prob- 
able model given the world's data. The RPR model has 
been shown to efficiently describe the p(j, K + )A observ- 
ables over a broad energy range [9[ HH [15] . The RPR 
framework, as it will be used in this work, has a mod- 
est number of adjustable parameters. The background in 
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the RPR framework consists of Reggeized K + (494) and 
* + (892) exchange in the t channel. With these assump- 
tions, the background part has two unknown phases and 
three unknown coupling constants. In addition to the 
background, the RPR model incorporates N* in the s 
channel. The improved version of the RPR model, as 
it will be introduced in this work, uses consistent N* 
interaction Lagrangians and this is an enormous asset 
in order to reduce the number of coupling strengths [H] . 
For each added N* one introduces one unknown coupling 
constant for J — \ and two unknown coupling constants 
forJ>§. 

The outline of the remainder of this paper is as fol- 
lows. In Section III Al the observables and kinematics of 
the N(fy,K)Y and N(e,e'K)Y reactions are introduced. 
In Section |IIB| we summarize the underlying assump- 
tions of the Regge-plus-resonance formalism used to de- 
scribe these reactions. Section |III| discusses a Bayesian 
approach to model selection. The computation of the 
Bayesian evidence is a high-dimensional problem which 
requires dedicated numerical methods and strategies. In 
Section |IIID| we provide details of these methods. A 
proof of principle of the adopted numerical strategy is 



described in Section HIE Bayesian methodology is ap- 
plied to determine the Reggeized background amplitude 
in Section |IV| In Section |Vj we determine the optimal 
resonant content for p(j, K + )A by evaluating a set of 
11 candidate resonances. For each of them we compute 
the relative resonance probability and the results are pre- 
sented in Section |V C| The results for the various pho- 
toproduction observables and predictions for electropro- 



duction observables are presented in Sections VD and 
|V E| A conclusion is given in Section [Vl| 



four-momenta are defined as follows: 

k* = (u*,k*) P* K = (E* K ,p* K ) 
p* N = (E* N ,-k*) p* Y = {E Y ,p* Y ) = {E* Y ,-p* K ). (2) 

The z-axis is the propagation direction of the incident 
photon, and the a;2;-plane is the reaction plane. 

Inserting the cm. momenta of Eq. ([2| into Eq. 
yields the expression for the unpolarized differential cross 
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where the transition amplitude can be written as the 

'A, 



product of the photon polarization vector e't and the 



hadronic current M\ nXy = £^ J XnXy . The hadronic 



current adopts the form 



JI nXy = u y Xy (Py)T» <(piv) 



(4) 



where t*L(pr) and u^ N (pjy) are the hyperon and nu- 
cleon spinors. 

The target (T) and recoil (P) asymmetries are defined 

as 
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da Xx =+h + da Xx =~^ 



(5) 



II. REGGE-PLUS-RESONANCE FORMALISM 
A. Observables and kinematics 

1. Photoproduction 

The unpolarized cross section for A ? "(7,if)F has the 
following expression 
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where v re x is the relative photon-nucleon energy, the 
p N (E N ,p N ), k(w,k), p K (E K ,p K ), and p y (E y ,Py) are 
the four-momenta of the nucleon, photon, kaon, and hy- 
peron. The A 7 , Ajv and Xy denote the photon, nucleon, 
and hyperon polarization. 

In the center-of-momentum (cm.) frame, the particles' 



where da = gfp^, and Ax is the nucleon and hyperon 
spin projection on the y-axis, respectively. The beam 
asymmetry S follows the definition 
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da — da" 



2da 
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where o" L (cr") refers to a linear photon polarization along 
the y(x) axis. 

Double polarization observables are defined as 



dcr(++) + da ( - 
da(++) + da(- 



-) + dai+~) + da(~+) ' 



(7) 



where (H — ) is a shorthand notation for (A^i = +sa, As = 
—sb), the polarizations of the particles A and B that 
determine the asymmetry. Beam-recoil pip/ , K + )A dou- 
ble polarization data are available for circularly (C x , C z ) 
and obliquely (O x ,O z ) polarized photon beams. These 
are more commonly expressed in the "primed" reference 
frame, which is rotated about the y = y'-axis over an an- 
gle Q* K , with 9* K the angle between the incoming photon 
and the outgoing kaon momentum in the cm. frame. 
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2. Electroproduction 



hadronic part, by defining the two tensors 



For incoming and outgoing electron four-momenta 
fei(ei,fci) and fofejfo) the electroproduction cross sec- 
tion in the one-photon exchange approximation (OPEA) 
has the following form 
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where the hadronic part of the reaction is evaluated in 
the 7* N cm. frame and the leptonic part in the labo- 
ratory frame (p^ b = (mjv,0)). This yields the following 
expression for the unpolarized differential cross section 
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where df 1 (e<j ) is the incoming (outgoing) electron en- 
ergy in the lab frame and W is the invariant energy. 

In the transition amplitude T\ , Ai and A2 are the 
polarizations (for high-energy electrons equal to the he- 
licities) of the incoming and outgoing electron. This am- 
plitude has a leptonic and a hadronic current, connected 
by a photon propagator 
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where the hadronic current J^ nXy is defined in Eq. Q, 
l^ 1 * 2 is the leptonic current, and k 2 = (k 2 — fci) 2 = — Q 2 . 

Therefore, T^\ 2 can be linked to M\ nXy and one can 
write 
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where the photon propagator was rewritten using the re- 
lation 



E (-^ 
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and Q 2 



-fc 2 is the photon virtuality. 

- Ai A2 



The tensor L^ lA2 is defined as a contraction between 
the photon polarization four-vector and the leptonic cur- 
rent 
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This allows one to replace the squared transition ampli- 
tude in Eq. ^ by 
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After this replacement, the separated cross sections or 
structure functions emerge. They do not depend on the 
kaon azimuthal angle <jf K and are defined as 
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where x - (16ff) 2 Wmjv r^fffa - 

Expressing Eq. ([9]) in terms of the separated cross sec- 
tions, we obtain 
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in which the dependence on <p* K has been made explicit. 
The virtual photon flux 
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and the virtual photon (transverse) polarization 



(22) 



(23) 



are defined in terms of the electron scattering angle e 
and the virtual photon three-momentum in the lab frame 

fclab- 



B. Regge-plus-resonance formalism 



Using Eq. ( 11 ) in the OPEA, one can conveniently sep- 



arate the quantum electrodynamics (QED) part from the 



This section deals with the dynamics of kaon produc- 
tion as described by the Regge-plus-resonance (RPR) 
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framework introduced in Refs. [H Oil EES] . The RPR 
model conjoins the economic description of high-energy 
data by means of Regge phenomenology with a single- 
channel hadrodynamical approach in the resonance re- 
gion. We have stressed the importance of the background 
diagrams in KY photoproduction for a correct determi- 
nation of the resonance parameters. In a isobar model, 
in which the amplitude is described as a sum of tree-level 
s, t and u-channel diagrams [HUH], the determination of 
the background is highly model dependent [20] . Another 
issue with isobar models is there violation of the Froissart 
bound [2"lTf2"5] . Indeed, the background amplitude of iso- 
bar models displays a power-law s a dependence at large 
energies where the exponent a depends linearly on the 
spin of the exchanged particles [3T]. The RPR approach 
overcomes these shortcomings by describing the nonres- 
onant contributions to the total amplitude by means of 
Regge theory [9]. 



1. Regge background 

Guidal, Laget and Vanderhaeghen showed that the ex- 
change of a limited number of Regge trajectories in the 
t-channel reproduces the high-energy, forward-angle data 
of both photoproduction [53] and electroproduction [2"5T - 
|2"T] of pions and kaons off the nucleon. Along those 
lines, the RPR background is obtained by Reggeizing 
the first materializations of the lightest kaon trajectories, 
K +(494) and K* +(892). The Reggeized amplitudes are 
obtained by replacing the t-channel Feynman propagator 
by a Regge one with the appropriate signature. 

The odd-spin and even-spin kaon trajectories are ob- 
served to coincide. The measured t-dependence of the 
da/dfl* K at large s does not display any pronounced 
structure, and this gives additional support to the strong 
degeneracy of the trajectories. Therefore, the Regge 
propagator reduces to 



^Rcggc( s ; *) 



' K 



sin (tvock (t)) \ e 
1 



T(l + a K (t)) 



(24) 



where a' K is the slope of the trajectory and the scale 
factor so is fixed at 1 GeV 2 . 

The relative sign between the odd-spin and even-spin 
propagators determines the phase (1 or e - i7rQ K(t)) f 
T^Rcrgc an d cannot be determined from first principles. 
The issue of determining this phase by comparing model 
predictions with data will be addressed in Section |IV| 

For vector mesons, we obtain the proper pole positions 
by subtracting the spin from the trajectories in the Regge 
propagator. For the K* + (892), the resulting propagator 



a K , (t)-ao 



i K , 



sin (-K(a K ,(t) - a )) 
■ ! *<"K" n f r(l + a(t)jf. -a ) 



(25) 



where a?o — 1. The employed parametrization for the 
#+(494) and K*+(892) trajectories is given by [20] 

a K (i) = 0.70 GeV" 2 (t - m 2 K ) , (26) 
a K . (f) = 1 + 0.85 GeV" 2 (t - m 2 K , ) . (27) 

2. Gauge restoration 

The K + (494) exchange diagram in the t-channel 
breaks gauge invariance. One way of restoring it is to 
add the electric part of the s-channel Born-diagram with 
the same coupling constant as in the K + exchange dia- 
gram. This procedure is also applicable for a Reggeized 
t-channel. It turns out to be essential for a proper de- 
scription of the forward-angle differential cross sections 
and of the beam asymmetries in charged pion photopro- 
duction [24] . 

For p(7, K )A, the gauge-restoring s-channel contri- 
bution is pivotal to account for the plateau in the differ- 
ential cross sections at very forward kaon angles or small 
\t\ [24]. The differential cross section for p(j, K + )A is 
shown in Fig. [T] 

Along similar lines, gauge invariance for electroproduc- 
tion can be restored by adopting the same Q 2 -dependence 
in both the electromagnetic coupling of the K + exchange 
diagram and the electric part of the s-channel Born term. 
In practice, this implies that a monopole kaon form fac- 
tor is assigned to the electric part of the proton exchange 
diagram |28j . This procedure has been shown to result 
in a reasonable prediction of the o\l/ot ratio [26 . 



3. Adding resonance contributions 

While Regge theory provides a fair description of me- 
son photoproduction observables at high energies and 
forward angles, there arc arguments that it can also 
be applied in the resonance region. Indeed, the notion 
of reggeon-resonance duality states that the amplitude 
should be reproduced by summing over all diagrams of a 
certain channel, be it the s, u or t-channel [29] . 

Even though the smooth s dependence of the Regge 
amplitude does not allow one to describe the structures 
in the resonance region, the global trends can be fairly 
reproduced |26j . Furthermore, the forward peaking of 
the differential cross sections supports large contributions 
from the nonresonant t-channel background. 

Inspired by these observations, Corthals et al. [H [El 
[T71 [TB] developed a hybrid model for KY photoproduc- 
tion dubbed Regge-plus-resonance (RPR). We will refer 
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to this model as RPR-2007. The RPR-2007 model uses 
amplitudes which consist of s-channel resonances and 
Reggeized ^-channel background terms. This approach 
has also been successfully applied to the electromagnetic 
production of tttt [30], as well as r] and 77' [31] . 

The Regge background amplitude of RPR-2007 is con- 
strained to above-resonance (y/s > 3 GeV), forward- 
angle (cos9* K > 0.35) data. By extrapolating the re- 
sulting amplitude to smaller y/s, one gets a parameter- 
free background for the resonance region. The s-channel 
resonances are coherently added to the background am- 
plitude, resulting in a hybrid amplitude for the resonance 
and high-s region. RPR-2007 describes the data for 
forward-angle photo- and electroproduction of K + A and 
K+TP [§1HH HE]. With regard to the 7V*'s, it includes 
the established PDG resonances Sn(1650), Pn(1710), 
Pi 3 (1720), the less established P 13 (1900), as well as the 
missing Di3(1900). The resonance parameters of the 
RPR-2007 model are constrained to the cos^ > 0.35 
data. 



0.5rr 




cos9 K 



Figure 1: (color online). The p(-y, K + )A differential cross sec- 
tion as a function of cos for the laboratory photon-energy 
bins cji a b = 1575 MeV, 1875 MeV and 2175 MeV. The line 
denotes the RPR-2007 result and the data are from refer- 
ences [32H35] . The RPR-2007 model is optimized against the 
cos,6* K > 0.35 data (indicated with the arrow). 



In Fig. [T]we confront the predictions of the RPR-2007 
model with a selection of differential cross-section data. 
At forward angles the data are nicely described, in stark 
contrast to the situation at backward kaon angles. This 
failure of the model at backward angles can be largely 
attributed to the adopted description for the spin-3/2 
resonance diagrams in RPR-2007 [TB]. Obviously, un- 
physical bumps at backward angles manifest themselves 
and the situation worsens with increasing lab photon en- 
ergy uJi a b. In the forthcoming Section it is pointed out 
how the introduction of consistent high-spin interactions 
can remedy this situation. 



4- Consistent high-spin interactions 

In the RPR-2007 framework, spin-3/2 resonances 
are described by the Rarita-Schwinger formalism [36] , 
Rarita-Schwinger fields, however, contain lower-spin 
components, which are not physical. In the nonintcr- 
acting Rarita-Schwinger theory these unphysical compo- 
nents are eliminated by imposing the so-called "Rarita- 
Schwinger constraints". These constraints, however, do 
not prevent the unphysical components from participat- 
ing in the interacting theory. The spurious lower-spin 
components generate non-localities, violate causality [37) . 
and must therefore be avoided. 

The spin-3/2 interaction Lagrangians that are used in 
the RPR-2007 model are inconsistent since they allow for 
the propagation of the unphysical spin-1/2 modes of the 
Rarita-Schwinger field. These Lagrangians involve the 
coupling of the spin-3/2 Rarita-Schwinger field through 
the so-called "off-shell tensor" , which contains a free pa- 
rameter. This off-shell parameter is associated with the 
unphysical contribution to the spin-3/2 interaction. The 
spin-3/2 resonance exchange diagrams of the RPR-2007 
model contain three off-shell parameters. 

In Ref. [TB] a consistent theory for the interaction of 
high-spin fermions was devised. There it was shown that 
an interaction theory that is invariant under the so-called 
"unconstrained Rarita-Schwinger gauge" is consequently 
a consistent theory, i.e. the unphysical components of the 
Rarita-Schwinger field decouple from a gauge-invariant 
interaction. 

In the updated version of the RPR framework, dubbed 
RPR-2011, the exchange of spin-3/2 resonances is de- 
scribed by the consistent interaction theory of Ref. [16] . 
In addition, the RPR model has been extended to include 
the exchange of spin-5/2 resonances. The expressions for 
the KYR(3/2) and KYR(5/2) interaction Lagrangians 
read [16] 

£kyr(3/2) = lfKYR 2 (3/2) %T^ Y d^ K + H.c, (28) 
m K 

£kyr( 5 /2) = - fKYI f /2) % y T'4, Y d^K + H.c. (29) 
% 

Here, ipY and 4>k represent the hyperon spinor and the 
kaon field. The factors fKYR{3/2) an d Ikyr(5/2) are 
strong coupling constants. Further, T = 1,T' = 75 for 
even parity resonances and T = 75, L' = 1 for odd parity 
resonances. The explicitly gauge-invariant fields and 
^ R v describe the consistent spin-3/2 and the spin-5/2 
resonances and read 

^ R = i(d^ u r R - (30) 

(31) 

where ip R an d ip^ denote the spin-3/2 and spin-5/2 
Rarita-Schwinger fields. The interaction Lagrangians for 
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the 7pi?(3/2) and jpR(5/2) couplings are given by 
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-(2) 

- 7J >R(3/2) 



^Za^rV^^+H.c, (32) 

(2) 
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H.c. (35) 



The electromagnetic tensor contains the photon field 
and is given by = d^A^ — d w A^. Further, ip p rep- 
resents the proton spinor and ^^3/2): K pi?(3/2^ K 



and 



(2) 
"pR(5/2) 



pfl(5/2) 

are electromagnetic coupling constants. 



The RPR-2007 model employs a Gaussian hadronic 
form factor (HFF) to regularize the transition amplitude 
beyond a certain energy scale. From the expressions (30) 



and (31 ) for the explicitly gauge-invariant fields, it is seen 



that the power of the momentum dependence of a consis- 
tent interaction rises with the spin of the exchanged par- 
ticle. In Ref. [Hj it is shown that unlike a Gaussian HFF 
a "multidipole-Gauss form factor" is capable of suppress- 
ing this momentum dependence. The functional form of 
this HFF reads 



F mG (s;m R ,A R ,T R ,J R ) = exp 



m 2 R ) 2 



A 4 



(Jr) 



Jr- 



m R f+m R Tl{J R ) i 



where T R (J R ) is defined as 



r R 

2^' 



(36) 



(37) 



In this expression, m Rl K R , T R , and J R denote the mass, 
the cut-off energy, the decay width, and the spin of the 
exchanged resonance. For J R = 1/2, Eq. (36) reduces to 
the familiar Gaussian HFF. The RPR-2011 model uses 



the multidipole-Gauss HFF of Eq. (36) in order to regu- 



larize the high-energy behavior of the consistent spin-3/2 
and spin-5/2 transition amplitudes. We use one common 
cut-off A R for all resonances. 



III. BAYESIAN INFERENCE 



Table I: Jeffreys' scale for the natural logarithms of evidence 
ratios A In Z = In [381 139] . It provides a translation be- 
tween the evidence ratio or Bayes factor and a qualitative as- 
sessment of the premise that model A is more probable than 
model B. 

\A\nZ\ < 1 Not worth more than a bare mention 
1 < \AhiZ\ < 2.5 Significant 
2.5 < I A In 2 1 <5 Strong to very strong 
5< lAlniJ Decisive 



A. Model comparison 

Using Bayes' theorem, P(A\B) P(B) = P(B\A)P{A), 
one can straightforwardly derive a quantity of interest 
for model comparison: the probability P(M\ {dk}) of a 
model M, given a set of experimental data {dk} 



P(M\ {4}) 



P({d k }\M)P(M) 
P({d k }) 



(38) 



The quantity P({dk] \M) is referred to as the marginal 
likelihood or the Bayesian evidence (Z). If the model 
M can have different outcomes, which are parametrized 
with a set of numbers CKjvr, marginalization yields 

Z = P({d k }\M) = J P({d k },a M \M)da M , (39) 
£{a M )ir(a M )da M . (40) 



This expression states that the Bayesian evidence is the 
integral of the product of two distributions: (i) the prob- 
ability of the dataset {d k }, given the set of parameters 
olm and the model M, and (ii) the probability of the set 
of parameters o.m, given the model M. The first fac- 
tor, P({d k } \ctM)M), can be identified as the likelihood 
function, C(cxm)- Any prior knowledge of the parame- 
ters' probability distribution before considering the data 
{dk} is contained in the second factor P(o.m\M), which 
is referred to as the prior distribution n(etM)- 

The quantity of interest for model comparison is the 
relative probability of a model M A versus a model M R , 
given the available experimental data {dk}. By applying 
Bayes' theorem (38), the evidence ratio or Bayes factor 



readily emerges from the expression for this probability 
ratio: 



P(M A \{d k }) 
P{M B \{d k }) 



P({d k } \M A ) P(M A ) 
P({d k } \M B ) P{M b ) 



iovP(M A ) = P(M B ). 



(41) 
(42) 



In this Section we outline how Bayesian inference can 
be used to constrain a framework like RPR against a set 
of data. 



The natural logarithm of the evidence ratio can be inter- 
preted qualitatively with the aid of Jeffreys' scale [3"5ll39| . 
given in Table [Tj 
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B. Probability of a resonance 

Bayesian inference can also be used to extract the phys- 
ical properties from the data. For example, does the fit 
to a set of photoproduction data provide evidence for 
the introduction of a hitherto unknown resonance? Wc 
present a procedure to calculate the relative probability 
of a certain nucleon resonance within a model for KY 
production, such as the RPR model. This procedure will 
help fill the need for an unbiased quantity that expresses 
the need for introducing an unknown resonance. 

Note that all probabilities mentioned in this subsection 
are implicitly conditional on a given framework. The de- 
pendence on the RPR framework AArpr is implied from 
now on, but will be omitted for the sake of clarity, i.e. 
P (X) = P (X | Mrpr). One can write the probability 
of a given resonance R, given experimental data {dk} as 



The quantity x ( a M) is defined as 



(43) 



Mi 



= J2 p ( R \ M h{dk})P(Mi\{d k }). (44) 

Mi 

The conditional probability P (R \Mi, {dk}) simply re- 
duces to one if the resonance R is included in the set of 
resonances Si used in the model variant M, , and zero oth- 
erwise. Therefore, the summation covers only a limited 
set of models, 

P(R\{d k })= £ P{M i \{d k }), (45) 

Mi\R£S z 

Mi\R€S t Vl JJ 

Applying Bayes' theorem, one finds that the evidence 
P({dk}\Mi) appears in equation (46). Assuming that 
there is no preference for any specific model before com- 
paring it to data, the factor pnyn is equal for all 
model variants i. Therefore, the factor can be omitted 
in all subsequent calculations for the probability ratios. 
This again reduces the calculation of relative probabili- 
ties P (Ri | {dk}) jP {R% | {dk}) to the evaluation of the 
evidence integrals of the form of Eq. 40 



C. Likelihood function 

Experimental data are usually reported to have nor- 
mally distributed errors and to be independent. The ad- 
dition of N squared normally distributed, independent 
random variables with mean and variance 1 results in 
a variable X = 

button mum 



XL=i x i that obeys a chi-square distri- 



x N ' 



2-1, 



-X/2 



2 w / 2 r(f ) 



(47) 



N IA 

y (dj 

8 = 1 



/i(ai)) 



(48) 



where N is the total number of data points, <7j is the error 
on data point di, and fiipiM) is the corresponding model 
prediction. The quantity x 2 ( a M) represents a sum of 
squares of normally distributed variables and is expected 
to obey the chi-square distribution of Eq. ( 47 1 . 

With a likelihood function of the form (47), we get the 
log-likelihood 



ln£(o!jvr) 



ln X 2 (a M )-^ln2 



lnT 



(49) 



where k is the number of degrees of freedom: this is equal 
to the number of data points N minus the number of 
free parameters. This correction is necessary because by 
constraining the free parameters using the data, one ef- 
fectively decreases the number of degrees of freedom. 

The x 2 { a M) and L{olm) are unknown functions of the 
model parameters clm and the numerical computation of 
the Bayesian evidence Z with the aid of the Eq. ( 40 ) in- 
volves a multidimensional integral J datM over the full 
parameter space. This is highly nontrivial from the nu- 
merical point of view. In the forthcoming section |III D| 
we outline the adopted strategy in order to compute the 
Bayesian evidence. 



D. Numerical computation of the Bayesian 
evidence 

For low-dimensional problems (d < 10), the Nested 
Sampling (NS) Monte Carlo algorithm by Skilling |4"2"ll4l] 
provides an efficient means to compute the Bayesian ev- 
idence. The posterior distribution P{olm\ {dk} , M) can 
also be computed by this algorithm. We employ this 
method to determine the Reggeized background ampli- 
tude of the RPR-2011 model @3]. The results of this 
analysis are reported in Section |IV| 

In high-dimensional problems, NS has been criticized 
for having a sharply decreasing acceptance rate as the 
likelihood constraint becomes more exclusive [32] ■ There- 
fore, high-dimensional problems call for an alternative 
numerical technique. If there is no need to determine the 
posterior distribution, or if the parameters are so-called 
"nuisance parameters", whose value are of no interest, 
other Monte-Carlo integration methods can be employed. 
One such method is the VEGAS algorithm by Lepage [3S] • 
VEGAS uses importance sampling: the points are sam- 
pled from a proposal distribution which approximates 
the normalized integrand. The proposal distribution is 
discretized in the form of an adaptive grid, of which each 
cell is sampled with an equal probability. This idea is 




Figure 2: (color online). The likelihood function C(a) (green 
surface) and the proposal distribution determined by VEGAS, 
represented by a grid (red lines) with a number of bins per 
dimension (here 8). 



illustrated in Fig. [2] The VEGAS algorithm is most suit- 
able if the integrand C{glm) can be approximated by a 
separable function. 

We have adapted the GNU Scientific Library (GSL) 
implementation [37] of the VEGAS algorithm to the in- 
tegrand of the evidence integral, which can assume very 
small values. This adapted VEGAS method, which we will 
refer to as log-VEGAS, requires a function which returns 
the natural logarithm of the integrand. The integral is 
computed while ensuring minimal loss of numerical accu- 
racy that would occur by exponentiation of this function. 
This measure is indispensable for the integration of small 
quantities such as a likelihood. 

Like any stochastic integration method, the log- VEGAS 
algorithm is apt to miss a highly localized maximum. 
We remedy this by locating the maximum with a genetic 
algorithm (GA) before performing the integration. Wc 
combine a rough search in the full parameter space using 
a GA and a subsequent fine search in a selected part of 
the parameter space using the gradient-based methods 
of MINUIT, the optimization module of the ROOT library 
[35] , This strategy has been successfully applied to a 
precise determination of resonance parameters by Ireland 
et al. (2]. 

The next step is to reduce the integration space to the 
volume around the peak, with a range of the order of 
three standard deviations in each dimension. The stan- 
dard deviation around the maximum can be calculated 
using the MINOS routine of minuit [35] , 

The first question that springs to mind is whether we 
do not risk underestimating the evidence by limiting the 
integration domain to the peak volume. We have ad- 



dressed this concern by applying this method to a toy 
example, which is detailed in the following Section. 



E. Toy example 



As a proof of principle, we apply the methods outlined 
in Section [lII D| to a tractable and realistic-sized problem. 
To this end, we use an event described by the function 
rrid(x) which is expressed in terms of a sum of d Legendre 
polynomials: 



m d (x) 



d-l 

E 

1=0 



aiPi{x) 



€[-1,1]. 



(50) 



The parameters ai(l — 0, .. . ,d — 1) are uniformly dis- 
tributed in in [—10, 10] and randomly generated. A mock 
data set with Gaussian noise is generated from rrid(x). 
The effectiveness of the GA is assessed by testing whether 
the values ai can be determined from the mock data. In 
the next step it is investigated whether Bayesian infer- 
ence can determine which model was used to generate 
a particular set of mock data. In essence, this amounts 
to use Bayesian inference to find the dimension d of the 
model from which the mock data are generated. 

We investigate the performance of a GA for models 
with a complexity ranging from d = 1 to d = 12. We 
consider 4000 data points, a size comparable to that of 
the world's K + A photoproduction data set. For each 
data set, we attempt to determine the parameters of the 
underlying model. We scale the population size in the 
GA linearly with d. Due to its random character, con- 
vergence times can vary greatly between the different GA 
runs. To account for this, we have repeated the GA 40 
times for each value of d, using a different, random set 
of parameters for each run. We have found that con- 
vergence occurs for all trial runs, and that the original 
parameters are reproduced by the GA with an error per 
parameter of the order of 0.5%. 

Can one determine the model which best describes a 
given data set from a number of model variants? This 
key question can be rigorously addressed using Bayesian 
inference. To illustrate the potential of this method, we 



generate mock data using the toy model of Eq. ( 50 1 at a 



fixed d. In a next step, we try to determine the under- 
lying model (including its dimension d) by calculating 
the Bayesian evidences for different trial models using 
the log-VEGAS method. This procedure is tested for data 
sets generated by models of different complexity: from 
d = 1 up to d — 12. 

The results of the log- VEGAS integrations over the en- 
tire parameter space are compared to those limited to 
the peak volume in Fig. [3j We find that the maximum 
evidence value corresponds to the correct model up to at 
least d — 12. 

Two striking conclusions can be drawn from Fig. [3j A 
first observation is that in the low-dimensional problems 
(d < 10), where the log-VEGAS result can be considered 
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Figure 3: (color online). The —\nZ values calculated with the log- VEGAS algorithm (red squares) and with a combined 
GA+MINUIT+Iog- VEGAS integration (blue diamonds), for different model dimensions. Each box corresponds with a model with 
dimension d, indicated with the grey band. The mock data set has 100 points. 



accurate, the two methods provide a comparable value 
for the computed integrals. This means that the like- 
lihood in the parameter space outside the peak region 
is small as compared to the maximum likelihood. Sec- 
ond, the results for high-dimensional models (d > 10) 
indicate that the bulk of the evidence is somehow over- 
looked by the global log- VEGAS integration. The global 
integrals for high-dimensional problems can be orders of 
magnitudes smaller than those that cover only the re- 
gion around the peak. This indicates that the search 
space for the log- VEGAS integration is too large in these 
high-dimensional problems, and a more dedicated search 
strategy is required. The results of Fig. [3] indicate that a 
combined GA+MINUIT+Iog-VEGAS integration strategy 
is apt to the task of dealing with high-dimensional prob- 
lems. 



IV. BACKGROUND SELECTION IN THE RPR 
MODEL 

In this section, we apply Bayesian inference to select 
the optimum model variant for the RPR background am- 
plitude. 



A. Parameters of the Reggeized background model 



The unknown phases in Eqs. (24) and (251 give rise to 
several model candidates. The possibility of the K + and 
K* + trajectories having a constant phase is excluded, 
as this gives rise to a recoil asymmetry P = 0, which 
disagrees with the data. In the forthcoming, the re- 
maining three possibilities, namely (rotating K + /rotat- 
ing K* + ), (rotating AT + /constant K* + ), and (constant 
K + /rotating K* + ), will be referrred as RR, RC, and CR. 



Apart from these variants, the background model has 
three continuous parameters proportional to the product 
of the strong and electromagnetic couplings, 



eG 



(51) 



Here, kk+k*+ is the transition magnetic moment for 
K* + (&92) -> 7A"+(494) decay. Further, the parameters 
feature the strong coupling constant gK+Ap °f the- K + 
trajectory and the tensor and vector couplings <?^*+ Ap 
of the K* + trajectory. 



Likelihood distribution and data 



As discussed in Section |III Cj the likelihood distribu- 
tion of the model parameters with regard to the data 
is the chi-square distribution of Eq. ( 47 ) . Recently, the 
CLAS collaboration published K+A [35 ] and K+E° @5] 
photoproduction data, featuring high-statistics differen- 
tial cross-sections and recoil polarizations. The data cov- 
ers nearly the full angular range and has 1.620 GeV 
< W < 2.840 GeV. The broad energy range makes it 
a great testing ground for both isobar, Regge and hybrid 
models such as RPR. Indeed, it includes measurements 
taken at energies up to W — 2.840 GeV, which is well 
above the resonance region. 

Sibirtsev et al. |11) demonstrated that the p(7,7T + )n 
and n(7,7r _ )p reactions display Regge-like behavior for 
invariant mass energies as low as 2.6 GeV. Further- 
more, Schumacher and Sargsian |50) pointed out that 
in the small- \t\ limit, the differential cross section for 
p(j, K + )A exhibits Regge-like scaling behavior cx s~ 2 
down to W « 2.3 GeV. One would therefore expect that 
a Regge background model optimized to the W > 3 GeV 
SLAC and DESY data [H], provides a fair description 
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of the W > 2.6 GeV CLAS data. However, this is not 
the case. Fig. [IJshows the W > 2.6 GeV K + A photopro- 
duction data, as well as the prediction of the Reggeized 
background model optimized to the W > 3 GeV data. 
Clearly, the Regge model overshoots the CLAS data by 
at least a factor of 2. There is an obvious discontinuity 
in the W dependence between the SLAC and CLAS data 
at cos 6* K w 0.865, 0.8, and 0.7. 

Dey et al. [52] showed that a small set of p{^,K + )K 
data from CEA [53] is inconsistent with the CLAS data. 
They find similar discrepancies between new CLAS data 
and old high-energy data from SLAC, DESY and CEA 
for other pseudoscalar meson production reactions. They 
conclude that there is a persistent normalization is- 
sue in the old high-energy differential cross-section data 
for a number of reactions, including p(j, K + )A and 
A" + )S°. The observations of Fig. Hadd support to 
these findings. 

Because of these observations, we opt to use the data 
from CLAS, which is consistent with other differen- 
tial cross-section measurements in the resonance region 
[531 [S3], in order to constrain the adjustable parame- 
ters in the Reggeized background model. We employ 
the statistical methods described in Ref. [44], using the 
2.6 GeV < W < 3 GeV CLAS data to compute the 
likelihood function. Below this energy region, resonance 
contributions become more important [11) . Because the 
validity of Regge theory is limited to small |t|, we use 
cos 9^ > 0.35 data to constrain the background parame- 
ters. With these criteria, we retain 132 differential cross 
sections and 130 recoil polarization P. This is over a fac- 
tor of four more data than for the combined SLAC/DESY 
data used in the analysis reported in Ref. [4*4"] . 



2. Prior distribution 

We opt to use a uniform prior distribution U for the 
coupling constants of Eq. (51). Under conditions of 
highly concentrated likelihood, compared to which the 



prior distribution varies mildly, the likelihood dominates 
the shape of the posterior distribution [32] • Accordingly, 
the evidence calculations will not be largely affected by 
the choice with regard to the shape of the prior distribu- 
tion. A sensitivity analysis will verify this assumption. 

The assumption that SU(3) symmetry is broken at the 
20% level yields the following prior ranges for gx+Ap [S3 



_ 45< 5k+Ap < 



47T 



-3.0 . 



(52) 



To our knowledge, for the K* + Ap vertex, no reliable the- 
oretical constraints are available [27]. We choose a uni- 
form distribution between —100 and +100 as the initial 
prior for (G v K , + ,G t K , + ) . To test the sensitivity of the 
results to the prior width, the calculations are repeated 
for a prior width of 2000 and 20000. 



3. Asymptotic behavior 

In the Regge (large s and small \t\) limit, one can ap- 
proximate s by — u for fixed values of t. This implies that 
the energy dependence of the cross section, which follows 
the power law s a ^' according to Regge theory, can be re- 
placed by (^) Q (*) [3D1 E7j. In an analysis of IV > 3 GeV 
data, this difference is not relevant, but at the energies 
considered here the difference between the two asymp- 
totic behaviors becomes noticeable. Therefore, we have 
investigated both options using Bayesian inference. 



B. Results 

1. Optimum background model variant 

The results of our analysis are listed in Table [Tl] The 
data clearly favor a model featuring an s Q< -*) dependence 
in the cross section and two rotating trajectories. In- 
deed, the difference in In Z with the second-best variant 
is 32.7 ± 1.4, which exceeds the value of 5 required for a 
decisive statement. 

The values of the coupling constants from the best 
model variant are 



9k+a p 



-3.6 ±0.3, 
9.0 ±0.5, 



V47T 

Gv 
K*+ 



G l K , + = 20.9 ± 0.4. 



(53) 



In comparison with the Bayesian analysis of Ref. [H] 
which was based on the W > 3 GeV data, the tensor 
coupling has changed sign, and its magnitude has 

decreased by about a factor of two. In Ref. [44. which 
was based on the W > 3 GeV data, the likelihood hyper- 
surface exhibited a distinct multimodal behavior, with 
different combinations of the coupling constants' rela- 
tive signs giving similar likelihoods. Interestingly, the 
increased amount of data used here causes the likelihood 
to be concentrated in only one quadrant of the param- 
eter space in (G V K , + , G f K , + ) and all sign issues for the 
coupling constants can be resolved. 

The expectation value ( 53 ) for gx+Ap is close to its 
SU(3) prediction of Eq. (52). Nevertheless, we have re- 
peated the analyses with a prior for gK+Ap broader than 
the condition of Eq. (52) in order to to test whether 



stronger SU(3) flavor symmetry breaking is compatible 
with the data. The results of this analysis are listed in Ta- 
ble [TTTJ By comparing the results of Tables [TT| and |III| one 
can conclude that the order of the models is unaffected by 
the operation of broadening the boundaries for the gx+Ap 
priors. Also the extracted value for gx+A P is n °t signifi- 
cantly affected by the broader limits on its prior distribu- 
tion. Its expectation value becomes g -^t^ p = —3.7 ± 0.3 
which is again compatible with the SU(3) value of —3.75. 
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Figure 4: (color online). The p(7, K + )A differential cross sections as a function of W for various cos 6^. The dashed lines 
represent the best model (RR) of Ref. [44 which follows from a Bayesian analysis of the W > 3 GeV data. The full lines 
correspond to the best model (RR) from Table [IT] optimized against the 2.6 GeV < W < 3 GeV CLAS data. The lines and 
the data are color coded according to cos 9^: from 0.4 (blue) to 1.0 (red). The orange lines correspond to cosf?^ = 0.95, the 
other lines have a value that corresponds to the CLAS cos#x bins, i. e. 0.865, 0.8, 0.7, 0.6, 0.5 and 0.4. Data are from Refs. 
[35] and IHI]- 



Table II: Logarithms of the evidence ratios ( A In Z = In (Z / Z max )) for the six model variants resulting from phase combinations 
and asymptotic behavior options in the two-trajectory Regge model for 7p — > K + A. The prior for the coupling constant g^+Ap 
is defined by the Eq. (521. The results are listed in order of decreasing probability for the lowest prior width, n = U{— 100, 100) 
for the G\ 



T K'+ 



couplings. 



K+ 1 K*+ phase 


asymp. 


TV = [/(-100, 100) 


tt = U(- 1000, 1000) 


7T = £/(-10000, 10000) 


RR 


s 


0.0 


0.0 


0.0 


RC 


s 


-32.7 ± 1.4 


-33.5 ± 2.7 


-31 ± 13 


RR 


(s - it)/2 


-359.7 ± 1.1 


-360.8 ± 7.2 


-389 ± 57 


RC 


(s - u)/2 


-432.9 ± 1.1 


-435 ± 8.9 


-472 ± 58 


CR 


s 


-2257.2 ± 1.1 


-2259.3 ± 6.4 


-2282 ± 31 


CR 


(s - u)/2 


-2425.5 ± 1.1 


-2426.3 ± 2.6 


-2440 ± 27 



2. High-energy predictions 

The high-energy differential cross section as calculated 
by the best model variant for p("f, K + )A is represented 
by the full lines in Fig. [4] As expected, the predictions 
are incompatible with the SLAC data. We attribute this 
to the normalization discrepancy discussed in Ref. |52j 
and Section [IV A 1 The polarization observables E [51] 
[55] and P [59. are not sensitive to normalization issues. 
Predictions for £ at = 16 GeV are shown in Fig. 5(a) 



for K + A. Fig. 5(b) shows the predictions for P at wi a b 



5 GeV. These predictions display an excellent agreement 
with data. By constraining the Reggeized background 
at 2.6 GeV < W < 3 GeV, one can predict P and S 
at W > 3 GeV. This highlights the predictive power of a 
Regge model at high W, and corroborates the assumption 
that the Reggeized background model can be constrained 
against W > 2.5 GeV observables. 

Summarizing the background evaluations, we find that 
the optimum two-trajectory Regge model for p(7, K + )A 
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Table III: Logarithms of the evidence ratios (A In Z = In {Z / Z max )) for the six model variants resulting from phase combinations 
and asymptotic behavior options in the two-trajectory Regge model for 7J3 — > K + A. A deviation of up to 40% from the 
SU(3) prediction for gx+A P i s allowed. The results are listed in order of decreasing probability for the lowest prior width, 
7T = t/(-100, 100) for the G^" + couplings. 



K + 1 K* + phase 


asymp. 


TV = U(-100, 100) 


tt = U(- 1000, 1000) 


7v = £/(-10000, 10000) 


RR 


s 


0.0 


0.0 


0.0 


RC 


s 


-17.8 ± 1.1 


-17.4 ± 2.9 


-15 ± 18 


RR 


(s - u)/2 


-359.7 ± 1.0 


-364.0 ± 16.0 


-387 ± 37 


RC 


(s - u)/2 


-432.7 ± 1.2 


-434.9 ± 5.3 


-474 ± 68 


CR 


s 


-2257.1 ± 1.2 


-2261.5 ± 7.7 


-2272 ± 28 


CR 


(s - u) /2 


-2425.6 ± 1.1 


-2426.3 ± 3.0 


-2426 ± 22 
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(a) 

Figure 5: (color online). Predictions of the best Regge model from Table (full red line) for the p(j, K + )A observables S and 
P at W > 3 GeV, as a function of —t. (a) E at cji ab = 16 GeV. Data are from Ref. 51. (b) P at LJi ab = 5 GeV. Data are from 
Ref. [59]. 



features two rotating phases and positive vector and ten- 
sor couplings. We also find that an asymptotic s a w 
dependence of the Regge amplitude is preferred over a 
((s — u)/2) a ^ one. This model will be referred to as 
Regge-2011 and determines the prior for the background 
amplitude of the RPR-2011 model. 



V. RESONANCE SELECTION IN THE RPR 
MODEL 

Given the world's p(j, K + )A data, this section ad- 
dresses the following questions: (a) From a proposed set 
of resonances, what subset features in the most probable 
model? (b) What is the probability of a proposed reso- 
nance Rl Bayesian inference allows one to answer these 
questions in a quantitative way. 



A. Data and resonances 

An overview of the available p(j, K + )A data is listed in 
Table IIVI In view of the normalization issue discussed in 
Section ITvl the data from SLAC ED and DESY GH are 



not included in the analysis presented below. The total 



number of data points which we incorporate is 6148, of 
which 3455 are differential cross sections, 2241 are single 
and 452 are double polarization results. We stress that 
after accounting for the error bars all data carry the same 
weight. 

We use the differential cross section data measured by 
the CLAS collaboration [HH [35] and LEPS [33]. Due to 
unresolved discrepancies with other data sets [BBJ, the 
SAPHIR differential cross-section data [32] is excluded. 
This decision is motivated by the fact that the different 
cross-section measurements by CLAS are internally con- 
sistent ;35] and consistent with the LEPS data [33] . To 
date, there is no independent measurement that confirms 
the SAPHIR data. 

The single polarization data consists of two sets of re- 
coil polarization data published by the CLAS collabora- 
tion [35] [B3], as well as a set from GRAAL [BT . The 
beam asymmetry data used in our analysis includes re- 
sults from LEPS (33] EH [60] and GRAAL [BT]. The in- 
cluded target asymmetries were determined by means of 
beam-recoil measurements by the GRAAL collaboration 
[63] . The included double polarization observables are 
beam-recoil asymmetries, consisting of C x and C z data 
by CLAS jBj)j and GRAAL 's measurements of O x > and 
O z , [B3J. 
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Table IV: Overview of the published experimental data for 
the reaction p(j, K + )A. 



Observable 


#data 


Experiment 


Year 


Reference 


da 


00 




iyoy 


Boyarski |51| 


79 n 


S APHTR 


9004 


niander R9l 




1 Q77 




ZUUD 


rsracuorci o4 




12 


LEPS 


2007 


Hicks [33] 




90fifi 

iUUU 




901 n 




E 


9 


SLAC 


1979 


Quinn 58 




45 


LEPS 


2003 


Zegers [60] 




04 


T FT"PQ 


zUUo 


Sumihama |54| 




4 




9007 


Hicks mi 




66 


GRAAL 


2007 


Lleres [61] 


T 


3 


BONN 


1978 


Althoff [52] 




66 


GRAAL 


2008 


Lleres [63] 


P 


7 


DESY 


1972 


Vogel [59] 




233 


CLAS 


2004 


McNabb [64] 




66 


GRAAL 


2007 


Lleres [61| 




1707 


CLAS 


2010 


McCracken [35] 


C x j Cz 


320 


CLAS 


2007 


Bradford 65 


o x , , o z , 


132 


GRAAL 


2008 


Lleres [63] 



tance of this resonance in the irN system was questioned 
in the most recent SAID analyses [68H70] . The Pu(1710) 
has also been identified in the ttttN system [71] . 

Furthermore, the two-star resonances Pi3(1900) and 
Fi 5 (2000)[87] are evaluated. The first of these, 
Pi3(1900), was found to couple to K + A by the Giessen 
group [75] and by the RPR-2007 model, and accounts for 
the structure in the energy dependence of the differen- 
tial cross-section data at W « 1900 MeV. Schumacher 
and Sargsian [50] show that the differential cross sec- 
tion data from CLAS 35 supports one or more reso- 
nances at W s» 2 GeV. Therefore, the consideration of 
the Fi 5 (2000) seems justified. The missing £>i 3 (1900) 
and Pn(1900) resonances earlier introduced in the Ghent 
isobar model QHHI], the RPR-2007 model 0E], and 
Kaon-MAID [1] are also evaluated. 

A conclusive statement with regard to the M ss 
1900 MeV resonances is extremely useful to improve our 
understanding of the nucleon's structure. Indeed, quark- 
diquark models do not predict a resonance at this energy 
[731174] . By contrast, a number of resonances with a mass 
of around 1900 MeV is predicted by CQMs [7S1[7B]. 



Table V: The nucleon resonances evaluated in the analysis 
given in the notation L2i,2j(M), along with their PDG status, 
spin (J) and parity (it), Breit-Wigner mass (M), width (F), 
and the uncertainty on the width (Ar). 



Resonance 


PDG status 


r 


M(MeV) 


r(MeV) 


Ar(MeV) 


511(1535) 




i/2- 


1535 


150 


±25 


Sn^O) 




1/2" 


1650 


150 


±20 


,015(1675) 




5/2- 


1675 


150 


-20/ + 15 


Pis (1680) 




5/2+ 


1685 


130 


±10 


£>i 3 (1700) 




3/2 _ 


1700 


100 


±50 


Pn(1710) 




1/2+ 


1710 


100 


-50/ ±150 


Pis(1720) 


**** 


3/2+ 


1720 


150 


-50/ ± 100 


#13(1900) 


missing 


3/2- 


1895 


200 




Pi 3 (1900) 




3/2+ 


1900 


500 


-360/ ± 80 


Pn(1900) 


missing 


1/2+ 


1895 


200 




Pis (2000) 




5/2+ 


2000 


140 


-40/ ± 30 



B. Likelihood function 

When calculating the likelihood function against a sin- 
gle data set, one usually does not take systematic errors 
cr sys into account. However, this course of action is not 
valid when multiple data sets are combined, as it would 
result in an underestimate of the likelihood. Assuming 
that systematic errors are independent and normally dis- 
tributed, the total errors can be determined by adding 
the systematic and statistical contributions in quadra- 
ture, 



(54) 



A more conservative estimate is to add the systematic 
and statistical errors linearly 



The 11 resonances considered in this work and their 
properties are listed in Table [V| We have "established" 
as well as "missing" nucleon resonances. As for their 
quantum numbers, mass, width and transition form fac- 
tors, we take the values quoted by the PDG. If these are 
not available, we employ the values determined by analy- 
ses based on CQM predictions [4]. This allows us to keep 
the number of adjustable parameters small. 

The established four-star resonances listed by the PDG 
are 5u(1650), P> 15 (1675), Pi 5 (1680) and P 13 (1720). 
The four-star Sii(1535) lies below the kaon production 
threshold, but is included because of its large decay width 
and its strong predicted coupling to the open strangeness 
sector [57] . To our knowledge, the contribution of the 
three-star Di 3 (1700) to p('f,K + )A is confirmed only by 
the Giessen analysis [5], The Pn(1710), which is found 
in some K + A analyses, is evaluated as well. The impor- 



CTstat ± <7 



sys 



(55) 



The numerical calculations for the Bayesian evidences are 
very demanding and it is prohibitive to run the calcula- 
tions with various choices for the values of Otot ■ In what 
follows, we outline an approximate method which allows 
one to relate the evidences computed with Eq. (54) to 
those which use Eq. ( 55 ) . 



Most often, the systematic errors a sys are computed by 
taking the squared sum of a number of partial systematic 
errors al ys from different sources. This approach is prone 
to underestimate the <7 sys . For a systematic error that is 
dominated by two errors with a comparable magnitude 
<7 sy s,i ~ o'sys^, one obtains in the conservative approach 



'sys 



"VM 



2cr. 



sys.l 



V2<7 : 



sys • 



(56) 
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The expression (54) for <r tot presupposes stringent inde- 
pendences between the various contributions. The a' tot 



of Eq. ( 55 ) provides a more conservative estimate of the 
evidence. In order to avoid overestimating the amount 
of information that is provided by the data, we will use 
a' tot in the forthcoming analyses. 



Figure 6: The integral over the likelihood function (black 
curve) can be approximated by the contribution for which 
C(a) « £ m ax, or conversely, for which x 2 ( a ) = Xmm (grey 
box). 



In a scenario where cr s 



cr S y S , the estimate (54 1 leads 



to <7t t ~ V^Cgtat an d to the following value for a more 



conservative estimate of a' tot 



°tot = ^tat + C. 



sys 



Cstat 



V2<T S 



1 + V 2 
V2 



"Ctot- 



(57) 



One can convert the Z values computed with the errors 



of Eq. (54) into a Z' which use a' tot . If the errors are 



multiplied by c, the log-chi-square distribution In C{olm) 
of Eq. ( 49 1 scales as 



In 



(k 



£'c(«Af) 
2)lnc- 



(58) 



Here, k denotes the number of degrees of freedom, Xr = 
X 2 /k is the reduced chi-squared as computed with the 
values tJtot ■ One can estimate the evidence resulting from 
the scaled likelihood function C'Jom) as follows. Insert- 



ing a uniform prior into the Eq. ( 40 ) yields the following 
expression for Z 



Z = 



C(a)da 



1 

A 



£(c*max)da, 



(59) 



where A = FJ. (Aaj) is the volume of the prior hy- 
percube. Indeed, if C(a) is the chi-square distribu- 
tion with xji( a ) f ar from its optimal value of 1 (e.g. 
x\ = 4), it falls rapidly with increasing x 2 («), and the 
bulk of the likelihood originates from a volume D where 
C(a) « £ max , or x 2 R [a) ~ X??„min> as illustrated in Fig.[6] 
By expressing the corrected likelihood C'Jot) in terms 
of the likelihood C(a) and the scaling factor of Eq. ( p38j ), 
the expression for Z' becomes 



Z' 



1 

A 
1 

A 

Z e 



£(«m. 



ix )da 

) p S(k, X 2 R (c 



c ) da 



-S , (fc,XB, mla .c) 



which yields our final result 

In Z' w In Z - S (k, x R ,min: c) 



(60) 
(61) 

(62) 



C. Identifying the resonance content of p(-y, K + )A 

The 11 proposed resonances of Table [V] give rise to 
2 11 = 2048 model variants. The Bayesian evidence Z' of 
Eq. (62) is computed for each model, resulting in a map 
of the RPR model space, shown in Fig. [7] 
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Figure 7: (color online). The evidences (—In Z') of the 2048 
model variants in the RPR model space (blue circles), as a 
function of the number of free TV* parameters. The latter 
value equals the number of A^* couplings plus one for the A_r. 
The smaller the value of — In Z' the higher the evidence. The 
best model for a fixed number of parameters is indicated with 
a red square. The model with the highest evidence, RPR- 
2011, is denoted with a black diamond. 



The parameters of the Reggeized background are as- 
signed localized priors of 20% around the values deter- 
mined in Section [TV) Therefore, the total number of ad- 
justable parameters is the sum of the number of N* cou- 
plings, of the three background parameters, and of the 
cut-off value A R of Eq. Q. The number of fitted N* 
couplings extends from 1 (one spin-1/2 coupling) to 18 
(4x1 spin-1/2, 4x2 spin-3/2, 3x2 spin-5/2 couplings). 
We adopt one common value for for all resonances 
with a uniform prior between 1.0 and 3.5 GeV. When se- 
lecting a prior distribution, it is good practice to ignore 
the data. Often an overestimation of the evidence re- 
sults from determining the likelihood and the prior with 
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Figure 8: (color online). The evidences (— \nZ') of the 2048 model variants in the RPR model space (blue circles). The purple 
diamonds correspond with the subset of models which contain the resonance indicated in the top right corner of each panel. 



a particular data set. In this work, the ranges of the 
prior distributions of the resonance couplings are selected 
on the basis of naturalness arguments. Indeed, the con- 
tribution of a single resonance is unlikely to exceed the 
total p(j,K + )A cross section (a w 5fib) by a large fac- 
tor. We performed calculations of the total cross sections 
(<7i{) in a model which includes a single resonance R and 
the Reggeized background. It is observed that the crite- 
rion aji < 25/i6 leads to absolute values of the coupling 
constants smaller than 100 in the adopted units conven- 
tion. Therefore, we adopt a uniform distribution for the 
priors of the resonance coupling constants in the range 
[-100,100]. 

Jeffreys' scale allows us to determine the "best" model 
from the 2048 variants. The model with the highest evi- 
dence has 14 TV* parameters (13 couplings and A/j) and 
features the 511(1535), 5u(1650), Fi 5 (1680), P l3 (1720), 



Pn(1900), Fi 5 (2000), and the missing D 13 (1900) and 
Pi3(1900). This model variant will be referred to as RPR- 
2011 [77]. The "second-best" model has two parameters 
less due to the absence of the Di3(1900). The difference 
in — In Z between the "best" and "second-best" models is 
2.3. This corresponds to significant to strong evidence in 
favor of RPR-2011. The difference with the other models 
is at least 6.8, which is consistent with decisive evidence 
for RPR-2011. 

In a next step, one can quantify the probability of 
each resonance separately by evaluating P(R\{dk}) of 
Eq. ( 46 ) . Fig. [8] visualizes which models are included in 
the sum. It is also instructive to calculate P (~ R | {dk}), 
the probability that a resonance is not required to de- 
scribe the reaction. The calculation of this quantity is 
completely analogous to Eq. (46). In Fig. [9| the mod- 
els which do not include a resonance R are visualized for 
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Figure 9: (color online). The evidences (— In Z') of the 2048 model variants in the RPR model space (blue circles). The orange 
diamonds correspond with the subset of models which do not contain the resonance indicated in the top right corner of each 
panel. 



each proposed resonance. 
The probability ratios 

In (P(R\{d k })/P(~R\{d k })), (63) 

are plotted in Fig. [10] A positive ratio indicates that 
the probability that the resonance R contributes to the 
reaction p^, K + )A is greater than the probability that 
it does not. Conversely, a negative ratio means that the 
data does not support the possibility that R contributes 
to the reaction. 

The results indicate that the resonances in RPR-2011 
are those that have a positive probability ratio. More- 
over, the two resonances with the highest probabilities 
are Pi3(1720) and Sii(1650). These are the two reso- 
nances that are also deemed important for the description 
of p(7, K + )A by most other models. There is decisive ev- 
idence that the P>i 5 (1675), Di 3 (1700), and Pn(1710) are 



not required to describe the p(-f, K + )A data. With re- 
gard to the resonance content in the 1800-2000 MeV mass 
range several suggestions have been in the literature, but 
no consensus has been reached. We have evaluated four 
states in that mass region (Table |v|): two have a two- 
star status and two are labeled as "missing". Our anal- 
ysis provides decisive evidence for three of these states: 
P 13 (1900), Pn(1900), and P 15 (2000). Note that the evi- 
dence for the "missing" Z?i3(1900) is significant to strong, 
but not decisive. 



D. Photoproduction with RPR-2011 

The p(7, K + )A observables presented in this subsec- 
tion are calculated with the RPR-2011 model parameters 
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fixed at their maximum likelihood values. The RPR-2011 
results are compared to the Reggeized background model 
Regge-2011 that was determined in Section |TV] From the 
difference between the Regge-2011 and the RPR-2011 re- 
sults one can infer conclusions about the role of the res- 
onances for the various observables. 

The p(j, K + )A differential cross section is displayed 
as a function of cos 0* K in Fig. 11 RPR-2011 provides 



a good description over a wide range of kinematics. For 
the lowest energies and backward angles there are devia- 
tions between the model and the data, hinting at possible 
missing dynamics such as u-channel contributions. It is 
striking that the i-channel background of the Regge-2011 
model already provides a reasonable description of the 
gross features of both the CJj & and cos 0* K dependence 
of the differential cross sections. The biggest effect from 
the resonance contributions is observed at the forward 
and backward kaon angles. 

Both RPR-2011 and Regge-2011 models exhibit a steep 
rise at extremely forward angles. At the three lowest 
u!i a b energies considered in Fig. [TTJ the inclusion of the 
resonances softens this rise and improve the goodness of 
the fit to the data. Note that the SAPHIR data (Fig. [1} 
suggest a plateau at forward kaon angles and that this 
feature is absent in the CLAS data. 

The angular dependence of the single-polarization ob- 
servables S, P, T is shown for three representative ener- 
gies in Fig. [12] The £, P, T receive stronger contributions 
from the A*'s than the differential cross sections. In con- 
trast to the high-energy situation considered in Fig. [5] the 
photon asymmetries in the resonance region are relatively 
small. The Regge-2011 reproduces the trend of increas- 
ing E with growing LOi a i,. The inclusion of the 7V*'s does 
not lead to a considerably improved quality of the fit. 
The recoil polarization P and target polarization T are 
highly sensitive to the resonance contributions. 

We now turn our attention to the double polarization 
observables. We stress that they represent but 7% of the 
total amount of data and that we give each data point an 
equal weight. As the bulk of the data is in the differential 
cross sections and to a lesser extent in the single polar- 
ization observables, the double polarization observables 
represent stringent tests of the RPR-2011 model. Per- 
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Figure 10: The relative resonance probabilities of Eq. |63| for 
each resonance listed in Table [V] 



haps somewhat surprisingly, Regge-2011 provides a good 
approximation to the double-polarization observables C x 
and C z . The observed trends C z ~ 1 and C x ~ C z — 1 [75] 
are well reproduced by both the Reggeized background 
(Regge-2011) and RPR-2011. This observation hints at 
the fact that the C x , C z are very background dominated. 
A large sensitivity to resonance contributions is observed 
for O x i and O x i , which are considerably better described 
by RPR-2011 than by Regge-2011. 



E. Electroproduction predictions with RPR-2011 

Electroproduction reactions have the potential to 
probe the hadron's electromagnetic substructure. The 
existence of this substructure can be parametrized by in- 
troducing a Q 2 dependence in the electromagnetic cou- 
pling constants. In addition, the p(e, e'K + )A reaction 
dynamics becomes sensitive to longitudinal couplings. In 
part, these couplings arise naturally from the photopro- 
duction amplitudes when gauge invariance is imposed. A 
peculiar class of longitudinal couplings vanishes for pho- 
toproduction reactions and cannot be constrained against 
real-photon data. The latter class of longitudinal cou- 
plings are neglected in the RPR model. In brief, we fix 
the basic reaction mechanism to the p(^, K + )K data and 
treat the electroproduction data as a test of the model. 
Thereby, we make some reasonable assumptions with re- 
gard to the electromagnetic form factors of the i-channel 
kaons and s-channel N* : s. Such an extrapolation of the 
p(j, K + )A amplitude to p(e,e' K + )A has been shown as 
reasonably successful for the RPR-2007 model [2]. We 
use the same N* helicity amplitudes (HA) as in Ref. [14] . 
Also the transition form factors for the spin-1/2 particles 
are those from Ref. [2]. The transition form factors 
for the spin-3/2 particles are derived from the consistent 
Lagrangians of Eqs. ( 32 ( — ( 33 ) . For the spin 3/2 and 5/2 



particles, the HA are calculated in the Bonn constituent 
quark model [76] . and the transition form factors are de- 
rived using Eqs. ( 34 )— ( 35 1. 

A comparison between recent low Q 2 = 0.030 — 
0.055 GeV 2 measurements and RPR-2011 predictions are 
contained in Ref. [53]. It was observed that RPR-2011 
provides a fair description of those data. Unseparated 
structure functions <tt + £&l at very forward kaon an- 
gles obtained in the 1970s are shown as a function of W 
and Q 2 in Fig. 14 together with Regge-2011 and RPR- 
2011 predictions. Obviously, at cos6>^- sa 1 the major im- 
pact of the intermediate resonances is to reduce the cross 
section by some modest factor. This is in line with the 
observations made for the real-photon differential cross 
sections of Fig. [TT] The electromagnetic form factors 
of the intermediate resonances reduce the effect of the 
A*'s with growing photon virtuality Q 2 . The RPR-2011 
model provides a fair prediction for both the Q 2 and W 
dependence of the data. 

Fig. [15] shows the energy dependence of the separated 
structure functions 07, and &t- In line with the data, 
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Figure 11: (color online). Angular dependence of the differential cross section at various incident photon energies uiab- The 
full red line represents the RPR-2011 model, the blue dashed line corresponds with Regge-2011. Data are from Refs. [32H35] . 



RPR-2011 predicts a <tl and ot of almost equal mag- 
nitude. The ut appears to be systematically underpre- 
dicted while ctl is somewhat overpredicted. The fair re- 
production of both the magnitude and the W dependence 
of ai, provides support for our assumptions with regard 
to the longitudinal couplings. 

Predictions for the transferred polarisation are pre- 
sented in Fig. [161 The Reggeized background model 



Regge-2011 as it was determined in Sect. IV B 1 predicts 
the flat W dependence and the magnitude of P' z ss 0.0 
and P' x R3 0.5. For P' z , the introduction of resonances 
worsens the quality of the agreement with the data ob- 
tained in Regge-2011. For P' x the effect of the N* is 
smaller and the quality of the agreement is better than 
for P', 



VI. CONCLUSION 

The RPR framework conjoins a Reggeized t-channel 
background with tree-level s-channel nucleon resonances 



from an isobar approach into an economical model for 
kaon photoproduction in and above the resonance region. 
The RPR model clearly separates nonresonant from reso- 
nant amplitudes which is an asset when searching for the 
properties of those (missing) resonances which contribute 
top(j 7 K+)A. 

We have used Bayesian inference to perform model se- 
lection both with regard to the resonant and nonreso- 
nant content of the RPR framework. It was shown that 
the Bayesian evidence Z is a quantitative measure for 
a model's fitness given data. The computation of Z 
requires involving multidimensional integrals which de- 
mand dedicated numerical methods. To that purpose 
we have proposed the "GA+MINUIT+Iog-VEGAS" inte- 
gration strategy. With this method one can reliably com- 
pute Z for models with a moderate number of adjustable 
parameters such as the RPR framework. 

First, the most probable model variant for the 
Reggeized background was determined against the 
2.6 GeV< W <3.0 GeV data. This involves the de- 
termination of three continuous and two discrete ad- 
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Figure 12: (color online). Angular dependence of the single polarization observables: beam asymmetry E (top row), recoil 
polarizati on P (middle row), and target asymmetry T (bottom row), at various incident photon energies ujuh- Line conventions 
as in Fig. [ll] Data are from Refs. [35l [60T [6TT [63l [64] . 



justable parameters. The extracted value for g^+Ap is 
compatible with the one predicted by SU(3) symme- 
try. Next, we have considered a set of 11 nucleon res- 
onances to determine the optimum resonant contribu- 
tion in the RPR p(j, K + )A framework. To this end, 
the Bayesian evidence was calculated for all 2048 model 
variants resulting from the various resonance combina- 
tions. The model with the highest evidence, dubbed 
RPR-2011, includes the resonances 511(1535), 511(1650), 
F 15 (1680), Pi 3 (1720), Pn(1900), F 15 (2000), Z>i 3 (1900), 
and Pi3(1900). An evaluation of the individual reso- 
nances' probabilities reveals that the two resonances with 
the highest evidence of contributing to p(j, K + )A are 
the Si i (1650) and Pi 3 (1720). There is decisive evidence 
that the Di 5 (1675), L>i 3 (1700), and Pu(1710) are not 
required to describe the current p(j, K + )A world's data. 
The computed evidence for the two-star Pi 3 (1900), the 
two-star Pi 5 (2000), and the "missing" P n (1900) is deci- 
sive, whereas for the "missing" £>i 3 (1900) it is significant, 
but not decisive. 

After fixing the basic reaction reaction mechanism to 



the p(j, K + )A data, the electroproduction data serve as 
a test of the model. In general our predictions for the 
electroproduction data are reasonably good which proves 
that the RPR-2011 model possesses predictive power 
and goes beyond a mere analysis framework. Therefore, 
we consider RPR-2011 as an efficient and robust model 
which can, for example, be used as an elementary pro- 
duction operator in strangeness production reactions in- 
volving the deuteron and finite nuclei. 
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Figure 13: (color online). Angular dependence of the beam-recoil double polarization observables with circular beam polar- 
ization, C x and C z (top rows), and with oblique beam polarization, O x i and O x i (bottom rows), at various incident photon 
energies wiab- Line conventions as in Fig. 11 Data are from Refs. [63) 165]. 
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